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Abstract. A theoretical framework for the determination of tidally induced radial- velocity variations in a compo- 
nent of a close binary is presented. Both the free and the forced oscillations of the component are treated as linear, 
isentropic perturbations of a spherically symmetric star. Resonances bet ween dynamic tides an d free oscillation 
modes are taken into account by means of the formalism developed by [Smeyers et al. (199^ ). The amplitude 
of the tidally induced radial-velocity variations seen by the observer depends on the orbital eccentricity and on 
the orbital inclination. The amplitude increases with increasing orbital eccentricity and is most sensitive to the 
value of the orbital inclination when 20° ^ i ^ 70°. In the case of a 5 Mq ZAMS star with a 1.4 M© compact 
companion, it is shown that resonant dynamic tides can lead to radial-velocity variations with amplitudes large 
enough to be detected in observations. The shape of the tidally induced radial-velocity curves varies from very 
irregular for orbital periods away from any resonances with free oscillation modes to sinusoidal for orbital periods 
close to a resonance with a free oscillation mode. Our investigation is concluded with an application to the slowly 
pulsating B star HD 177863 showing the possibility of resonant excitation of a high-order second-degree g^-mode 
in this star. 
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>- 1. Introduction 
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In close binary systems of stars, each component is sub- 
ject to the time-dependent tidal force exerted by its com- 
panion. A commonly used approach for the study of tidal 
effects in close binaries is based on the expansion of the 
tide-generating potential in Fourier series in terms of mul- 
tiples of the companion's mean motion. Through these 
expansions, the tidal action of the companion induces an 
infinite number of partial dynamic tides in the star, each 
with its own forcing angular frequency. 

In binaries with shorter orbital periods, the forcing 
angular frequencies may be close to the eigenfrequen- 
cies o f the free osc illation modes g'^ of the component 
stars ( [Cowling 1941 ) . These proximities lead to resonances 
which enhance the tidal motions of the mass elements and 
can have significant consequences for the observational 
properties of the binary. 



Rocca (1987) to include the effects of the Coriolis force in 
slowly rotating binary components. 

The luminosity variations associated with tidally ex- 
cited oscillation modes in close binaries were investigated 
by Kumar et al. (1995| ). The authors proposed a new 



method for the determination of the orbital inclination 
based on fitting the shape of theoretically derived light 
curves to the shape of observationally determined light 
curves. They applied their results to various polytropic 
stellar models and to the binary pulsars PSR J0045-7319 
and PSRB1259-63. 



Zahn (1970) studied resonances of dynamic tides with 
low-frequency g^-modes in non-rotating stars consisting 
of a convective core and a radiative envelope. For his pur- 
pose, he used an asymptotic representation for the eigen- 
functions of the g^-modes which he established neglecting 
the perturbation of the gravitational field due to the star's 
tidal distortion. Zahn's treatment was later generalised by 



More recently, smeyers et al. (1998) (hereafter referred 
to as Paper I) derived semi-analytical solutions for the 
components of a resonant dynamic tide by means of a 
two-time variable expansion procedure in which both the 
free and the forced oscillations of the star were treated as 
linear, isentropic oscillations of a spherically symmetric 
star. The authors concluded that, at the lowest-order of 
approximation, the oscillation mode involved in the res- 
onance is excited with the forcing frequency of the reso- 
nant dynami c tide. Their investiga tion was subsequently 
extended by Willems et al. (2001 ) to take into account 
the effects of radiative damping in the nonadiabatic sur- 
face layers of a star. 
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The duration of a resonance is determined by the 
combined effect of stellar and orbital evolution. Witte & 



Savoni e (1999, 2001) have shown that when both effects 
are taken into account, a dynamic tide can easily become 
locked in a resonance for a prolonged period of time. Such 
a long-term resonance can have a significant influence on 
the secular evolution of the orbital elements and increases 
the probability of detecting the tidally induced oscillations 
in the star. 

From an observational point of view, firm evidence 
of tidally excited oscillations is still scarce, possibly due 
to the absence of any systematic observational study on 
the nature of the pulsations discovered in close binaries. 
Studies of i ndividual st ars are repo rted in the literature 
e.g. crSco ([Fitch 1967|), 14 AurA ( |Fitch fc Wisniewski 



1979|)7|Spica ([Smith 198"5a|,|h|), 1 6 Lac (|ChapeUier et al 



1995}^V539 Arae ( Clausen 1996 ), among others, but the 
resonant nature of the modes was hardly ever proven. An 
introduction to a systematic observational study in the 
case of early- type binaries is given by Aerts et al. (199^ ) . 
Their initiative resulted from a more general long-term 
systematic study of line-profile variations in /? Cep stars 
and in slowly pulsating B stars, which led to the assess- 
ment that many of these types of pulsators belong to 
a close binary. Further evidence for the presence of B- 
type y-mode pulsators in close binaries was found by De 
Cat et al. (200C| ). An additional systematic observational 



project in this respect is the SEarch for FOrced Nonradial 



Oscillations (SEFONO) by [Harmanec et al. (1997[ ). 

In this investigation, our aim is to determine the radial- 
velocity variations associated with resonantly excited os- 
cillation modes in binaries with short orbital periods. To 
this end, we use semi-analytical solutions for the compo- 



nents of a resonant dynamic tide derived by Smeyers et al, 
(1998|)] We conclude our investigation by applying our re- 
sults to the slowly pulsating B star HD 177863. 

The plan of the paper is as follows. In Sects. 2 and 3, 
we present the basic assumptions adopted in our investi- 
gation and we decompose the tide-generating potential in 
terms of spherical harmonics and in Fourier series in terms 
of the companion's mean motion. In Sects. 4 and 5, we de- 
rive an expression for the radial- velocity variations asso- 
ciated with a resonant dynamic tide. The contributions of 
non-resonant dynamic tides to the tidally induced radial- 
velocity variations are determined in Sect. 6. In Sect. 7, 
we determine the total variation of a star's radial veloc- 
ity due to the various resonant and non-resonant dynamic 
tides and we illustrate the role of the orbital eccentric- 
ity and the orbital inclination. In Sect. 8, we apply our 
results to a binary consisting of a 5 Mq zero-age main se- 
quence star and a IAMq compact companion. In Sect. 9, 
the possibility of a resonantly excited oscillation mode in 
the slowly pulsating B star HD 177863 is investigated. The 
final section is devoted to concluding remarks. 



2. Basic assumptions 

Consider a close binary system of stars that are orbiting 
around each other in an unvarying Keplerian orbit with 
semi-major axis o and orbital eccentricity e. The first star, 
with mass Mi and radius i?i , is rotating uniformly around 
an axis perpendicular to the orbital plane in the sense of 
the orbital motion. The angular velocity is assumed to 
be low so that the effects of the Coriolis force and the 
centrifugal force can be neglected. The second star, with 
mass M2, is treated as a point mass. 

We start from the right-handed orthogonal frame of 
reference Cix'^x'^x'^ introduced in Paper I. The origin of 
the frame of reference coincides with the mass centre Ci 
of the uniformly rotating star and the x'^x'^-plane corre- 
sponds to the orbital plane of the binary. The directions 
of the x'^- and the a;''^-axis coincide with the direction 
from the star's mass centre to the periastron in the com- 
panion's relative orbit and to the direction of the star's 
angular velocity fl, respectively. 

As in Paper I, we pass on to an orthogonal frame of ref- 
erence whose origin and a;"^-axis correspond to the origin 
and the x'^-axis of the frame of reference Cix'^x'^x'^ , but 
whose x"^-axis and x"^-axis are corotating with the star. 
We use this frame of reference to describe the nonradial 
oscillations induced in the star by the orbiting compan- 
ion. With respect to the corotating frame of reference, 
we introduce the spherical coordinates r = (r,9,(j)). The 
transformation formulae from the Cartesian coordinates 
x'^, x'^, x'^ to the spherical coordinates r, 9, cj) are given 

by 



x'^ = r sin 9 cos ( 
x'^ — r sin 9 sin ( 
x'^ — r cos 9. 



(1) 



In order to study the radial-velocity variations due to 
the tidal motions of the mass elements located at the star's 
surface, it is convenient to pass on to a frame of reference 
whose origin coincides with the mass centre of the tidally 
distorted star and whose polar axis is directed from the 
mass centre to the observer. We perform this transforma- 
tion in two steps. 

First, we introduce an orthogonal frame of reference 
Cio^a^a^ whose origin and a'^-axis coincide with the ori- 
gin and the cc'^^-axis of the frame of reference Cix'^x'^x'^ . 
The direction of the a^-axis corresponds to the direction 
from the star's mass centre to the ascending node in the 
relative orbit of the companion. 

Next, we pass on to an orthogonal frame of reference 
Cia'^a'^a'^ with the same origin and the same a'-'^-axis 
as the frame of reference Cia^a^a"^ , but with its a''^- 
axis directed from the star's mass centre to the observer. 
The Cartesian coordinates a'^, a'^, a'^ are related to the 
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written as 




observer 



Fig. 1. The geometry involved in the introduction of the 
various frames of reference. The ascending node and the 
periastron of the relative orbit of the companion are de- 
noted by the points labeled A and P, respectively. 
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(2) 



Here to is the longitude of the periastron in the relative 
orbit of the companion, and i is the inclination of the 
orbital plane with respect to the plane perpendicular to 



the line of sight (see, e.g., Green 1985). With respect to the 
frame of reference Cia'^a''^a''^ ^ we introduce the spherical 
coordinates r = {r' ,0' , (f>'). 

The geometry involved in the introduction of the var- 
ious frames of reference is illustrated in Fig. ||. 



St W (r, t) = St 
exp [i ((Tt t — hnr)] 



EE 

£=2 m=-e k=~oo 



Wl^m,k {r ) 



(4) 



where ctt = kn + mCl is the forcing angular frequency 
with respect to the corotating frame of reference, r is a 
time of periastron passage, and 



GMi 



r 



YT 



(5) 



Here G is the Newtonian constant of gravitation and ce^rn,k 
a Fourier coefficient determined as 

Ce,m,k = 



1 



ie-\m\) 

{i+\m\) 
(1 + e cos 



pi™i(o) 



1 



^ cos(fc M + mv) dv. 



(6) 



In the latter expression, Pp{x) is an associated Legendre 
polynomial of the first kind, and M and v are respectively 
the mean and the true anomaly of the companion in its 
relative orbit. 

The Fourier coefficients C£^m,k decrease with increasing 
values of the multiple k of the companion's mean motion. 
The decrease is slower for higher orbital eccentricities, so 
that the number of terms that has to be taken into account 
in Expansion of the tide-generating potential increases 
with increasing values of the orbital eccentricity. For i = 2, 
the only non-zero coefficients C2,m,fe are those associated 
with m = —2,0,2. A more elaborate discussion of these 
coefficients can be found in [Willems (2000| ). 

For realistic orbital periods, the forcing angular fre- 
quencies ctt may be close to the eigenfrequencies of the 
star's free oscillation modes g^. The tidal action exerted 
by the companion is then enhanced and the oscillation 
mode involved is resonantly excited with the forcing fre- 
quency of the dynamic tide (Paper I) . The effects of reso- 
nant excitations of modes in close binaries are particularly 
important in the cases of resonances of dynamic tides with 
(7+-modes of a lower radial order. 



3. The tide-generating potential 

The tidal force exerted by the companion is derived from 
the tide-generating potential eTW{r,t), where Et is a 
small dimensionless parameter defined as 



(3) 



RA M2 
a ) Ml ■ 

This parameter corresponds to the ratio of the tidal force 
to the gravity at the star's equator. 

Following Polflict fc Smeyers (1990 ), we expand the 
tide-generating potential in terms of unnormalised spher- 
ical harmonics Y"^{6,(/)) and in Fourier series in terms of 
multiples of the companion's mean motion n = 27r/Toi.b, 
where To^b is the orbital period. The expansion can be 



4. Tidal displacement due to resonant dynamic 
tides 

We consider a single partial dynamic tide generated by 
the term 

£TWe,m,k{r)exp[i{aTt- knr)] (7) 

in Expansion (||) of the tide-generating potential and as- 
sume the forcing angular frequency ax to be close to the 
eigenfrequency ae,N of a free isentropic oscillation mode 
of radial order N that is associated with the spherical har- 
monic ¥^{9,4)). For the sake of simplification, we intro- 
duce the abbreviation S = {£, m, N) for the wave numbers 
£, TO, and N identifying the oscillation mode involved in 
the resonance. 
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Furthermore, let the relative frequency difference 

(8) 



In addition, the spherical harmonics 



are 



e = 



be of the order of the ratio ex of the tidal force to the grav- 
ity at the star's equator. At the lowest-order of approxi- 
mation, the components of the resonant dynamic tide with 
respect to the local coordinate basis d/dr, d/d6, d/dcj) are 
then given by 



iSr)T (r,t) 



1 St 
1 ~2 



P\ (cos 9) exp [\ [m (f) + (Tt t — k n t)] . 

{00)t [r, t) = - — Ci^rn,k Qi,N 5 

e 2 r 



9Pj"'(cos( 

m 

iS^)T (r,t) = 



— exp [i (m (p + (Tt t — k n t) 

TT Ce,m,k Wl,N 5 



(9) 



. pl™l (cos 9) exp [{{mcj) + ~ krir)] 

sin 9 

[Paper I, Eqs. (57)]. Here $£,Ar(r) and r]e^]y{r) are the ra- 
dial parts of the radial and the transverse component of 
the Lagrangian displacement of the oscillation mode S 
[Paper I, Eqs. (43)], and Qi^n is related to the work done 
by the tidal force through the oscillation mode S [Paper I, 
Eq. (50)]. 

The components of the resonant dynamic tide with re- 
spect to the local coordinate basis d/dr', d/d9', d/dcf)' are 
obtained from Solutions (^ by application of the trans- 
formation formulae for contravariant vector components. 
The transformation formula takes the form 



where 

dr' _ ^ 

dr ' 
89' 

'89 
89' 



( 8r' dr' 8r' \ 

dr 89 dcf) 

89' 89' 89' 

Ih ~d9 ^ 

8(j)' 8(f>' 8(1)' 

Va7 '89 '86 J 



( {6r)^\ 
\{54>)t) 



(10) 



= sm I cos 



D{9',cl)') 



sm I cos ( 



89 



sin^ 9'D{9',(t>'y 
cos i — sin i sin <b' cot ( 



(11) 



and the other partial derivatives are equal to zero. The 
function D {9' , (f>') is given by 



D 



cos (f)' + (cosi sm(f)' — s'mi cot 9'y 



1/2 



(12) 



transformed according to 
Yf^{9,(t)) = exp[-im(cj + f7t + 7r/2)] 



^ a,,™^fc.(z)r/(0', (/.'), (13) 

k' = -i 

where the coefficients a£,m,fc'(*) functions of the orbital 
inclination defined as 



ae,m,k'ii) = Mm^ ^^k' ie + my.{i~ k')\ exp(ifc'7r/2) 



E 



{-ly-^-l (cOS: 



* 27'+m+fc' / . \ 2l — 2j — m — k 

sniTT 



j\ {£ - m - j)! {e-k' - j)! (m + fc' + j)! 



(14) 



(Jeffreys 1965). In this definition, the summation is per- 
formed over all values of j satisfying the relations 

j ^ 0, j > — m — k' , j < £ — m, j < £ — k' , 

and the factors /x^' are given by 



IJ-k' = 1 



for k' > 0, 
< 0. 



(15) 



At the lowest order of approximation, the components 
of the resonant dynamic tide with respect to the local 
coordinate basis d/8r', 8/89', 8/dcj)' then take the form 



1 St 



{Sr')T (r', t) = - — Ct,ra,k Qi,N ^i,N{r') 

f 

\k'\ 



^ at,m.,k'{i)Pi '(cose') 



k'=-e 
exp 



\^i\^k' (j)' + kM -m(^Lj + ^^ j, 
{69 )t (r ,t) ^ - — Ci^rn,k Qi,N ;5 — 



ai,m,k'{'i) 



8Pf\cos9') 
89' 



k'=-e 
exp 



[i[k'<p' 



kM - miuo + 



{0<P )t [r ,t) ^ - — Ci^m,k Ql,N ,2 ■ 2ai 

s 2 r'^ sin'^9' 



a£,m,fc' {i) i k' Pj'' ' (cos 9') 

li^/cV' + fcM-TO^LJ+l) I 



k'=-e 
exp 



(16) 



5. Radial-velocity variations associated with 
resonant dynamic tides 

The tidal motions of the mass elements located at the 
star's surface contribute to the radial-velocity variations 
seen by the observer. Since the mass elements of the undis- 
torted equilibrium star are assumed to be at rest with 
respect to the corotating frame of reference, the radial 
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component of the velocity field associated with the tidal 
motions is given by 



vt 
where 



dt 



(-Cq/s) = iar ((5r)j, • (-6^/3) 



e„'3 = cos ( 



, d sin 6*' d 



(17) 
(18) 



dr' r' 80' 
is a unit vector in the direction of the a''^-axis. 

In the case of a resonance between a dynamic tide and 
a free oscillation mode, the tidally induced radial- velocity 
variations due to the mass elements located at the star's 
surface take the form 

1 Et 



e 2 



k'=-e 



cose' ^e^N{Ri) '(cos6'') 



Ri 



89' 



exp \ T k (p + k M — m (u! + 



The velocity variations seen by the observer are ob- 
tained by taking the average of Eq. ( |l9| ) over the visible 
hemisphere of the star. When one neglects the distortion 
of the stellar surface and the perturbation of the specific 
intensity due to the tidal motions of the mass elements, 
the average is determined as 

VT;res{t) = 



/o^" lo^^ ^ (^') WT;rcs(i?i, f', <t)']t) sin0' COS 9' de'd(j)' 



(20) 



!T lo^^ ^ (^') sin ^' cos e'd9'd(j)' 
where h {9') is a limb-darkening law of the form 
h{9')^l-u + ucose', ue[0, 1]. (21) 

After performing the integration over the azimuthal angle 
6' , one obtains 



1 Et 

7 T 



i]e,N{Ri) 



Ri 



Q 



exp < 1 



kM -m[Lj + 



m 



(22) 



Here f^j(u) and frj/{u) are functions of the limb- 
darkening coefficient u defined as 



/ h {9' ) Pt^ {cos 9') sin 9' cos"^ 9' d9', 
Jo 



7r/2 



^"^"'^'^ sinH' cos9'd9'. 



89' 



(23) 



(24) 



The functions are described in more detail in Appendix 

In addition to the partial dynamic tide with forcing an- 
gular frequency cr^, we also take into account the partial 



dynamic tide with forcing angular frequency — or, which 
is also resonant. By the use of the symmetry properties 
Ci-rn-k = Q,m,fc and ai-rnM = af,m,o, the global solu- 
tion for the variation of the star's radial velocity due to a 
resonance of a dynamic tide with a free oscillation mode 
takes the form 



VT;rcs(i) = £T af,m,o(«) Q,m,fc(e) Vf.,.^^{t), 

where 

1 



(25) 



■ (Tt 



Vi.NiRi) 



h,N sm 



kn{t t) — m (^Lo + — ^ 



(26) 

In this expression, the products Qe,N £,e,N{Ri) and 
Qe.N ''le,N{Ri) are independent of the normalisation 
adopted for the determination of the free oscillation mode 
of degree £ and radial order N that is involved in the res- 
onance. They also decrease rapidly with increasing values 
of N so that the influence of higher-order modes on the 
tidally induced radial-velocity variations is substantially 



smaller than that of lower-order modes (Willems et al 



Ui^j 1997|; [Willems 2000|). 



6. Radial-velocity variations associated with 
non-resonant dynamic tides 

We again consider a single partial dynamic tide generated 
by a term of the form 

£TWe^rn,k{r)exp[ii<7Tt- knr)] (27) 

in Expansion (^ of the tide-generating potential, but now 
assume the forcing angular frequency (Tt to be far from 
any of the eigenfrequencies associated with the star's free 
isentropic oscillation modes. The components of the tidal 
displacement field generated by this term can be expressed 
as 



{Sr)T (r, t) = St Ci^m,k ^rir) 

pj™' (cos 9) exp [i (m (j) + ctt t — k n t) 



{S9)t ir,t) = £TC£,,„,fc 



9pl"i(cos( 



89 

{S<p)T{r,t) ^ ETCl^rn.k ^2 
,|m| 



exp [i (m (j) + axt — knr)] 



' (cos 9) exp [i (m ( 



sin2 6' 
(Tt t — fc n r)] 



(28) 



where S.Ti'f) and rixir) are solutions of the homogeneous 
fourth-order system of linear differential equations 



dr 6 

drjT 
dr 
1 d 

j,2 dj. 



1 



d^j 
dr 



4t H 

9 



Vt - 



VT 



(7. 



1 

— Wi 



9 



*T,(29) 



(30) 
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(31) 



Here ^'t(?') is the radial part of the total perturbation 
of the gravitational potential (for a definition see, e.g., 
Polfliet fc Smeyers 1990| ), g the local gravity, the square 



of the isentropic sound velocity, and the square of the 
Brunt- Vaisala frequency. The functions ^t(''), VTir), and 
^xir) depend on the values of the azimuthal number m 
and on the Fourier index k through the forcing angular 
frequency ctt- 

The solutions of Eqs. (|2|) - (|3l|) must satisfy bound- 
ary conditions at the star's centre and at the star's surface. 
At r = 0, the radial component of the tidal displacement 
must remain finite. At r = the Lagrangian perturba- 
tion of the pressure must vanish, and the continuity of the 
gravitational potential and its gradient requires that 



dr 



= -{2£+l) 



£ + 1 
Rl 

GMi 



Rl 



(32) 



By comparison of Eqs. (^8|) with Eqs. (|9|), it follows 
that the contribution of a non-resonant dynamic tide to 
the tidally induced radial-velocity variations takes the 
form 

VT;e,m,kit) = et ai,mfi{i) ci^m,k{e) Vf.^„^f.{t) (33) 
with 

^T;£,m,fc(^) = 

1lT:e,m,k{Rl) 



2 Ctt /e,£(u)CT;£,m,fe(^l) " /^,^(w) 

kn{t — t) — m + — ^ 



Rl 



(34) 



7. Total variation of a star's radial velocity due to 
its tidal response 

In accordance with Expansion ^ of the tide-generating 
potential, the total variation of the star's radial velocity 
due to the tidal motions of the mass elements located at 
its surface, is obtained by adding the contributions stem- 
ming from the various resonant and non-resonant partial 
dynamic tides. 

When one restricts the contributions to those resulting 
from the tides associated with the second-degree spherical 
harmonics, the total variation of the star's radial velocity 
due to its tidal response can be written as 



oo 

^T;rad(0 = ^ 1 02,0,0 («) C2,0,fc(e) V^T;2,0,fc(0 



(35) 



k=l 



+a2,2,o(j) [c2,-2,fc(e)T^T;2 -2,fe(0 + C2,2,fc (e) V^y.a^a.fe )] } . 
where 



a2,o,o(i) = 2 « - l) ' 

a2,2,o(«) = a2,-2,o(«) = 3 siv? i. 



(36) 



In the case of a resonance of a dynamic tide characterised 
by the azimuthal number m and the Fourier index k 
with a free oscillation mode of radial order iV, the fac- 
tor „ k{t) must be replaced by the factor Vf.^^^{t). 

The contributions of the various partial dynamic tides 
to the observed tidally induced radial-velocity variations 
depend on the orbital inclination and on the orbital eccen- 
tricity through the coefficients a2,m.o(i) a-nd C2,m,fc(e), re- 
spectively. The variation of the product of the coefficients 
0'2,mfi{i) and C2,m.fc(e) is displayed in Fig. ^ as a function 
of fc, for the orbital inclinations « = 5° and i = 85°, and 
the orbital eccentricities e = 0.3 and e — 0.5. The product 
decreases with increasing values of k due to the rapid de- 
crease of the Fourier coefficients C2.m,fc(e) for higher-order 
harmonics in Expansion (^ of the tide-generating poten- 
tial. 

In the case of the lower orbital inclination i = 5°, 
the products of the coefficients a2.m,o(*) and C2.m,fe(e) are 
largest in absolute value for the coefficients associated 
with the azimuthal number m — Q. The prevalence is more 
pronounced for the eccentricity e = 0.5 than for the ec- 
centricity e — 0.3 due to the increase of the coefficients 
C2,o,fc(e) with increasing values of the orbital eccentricity. 
For larger values of the orbital inclination, the products 
of the coefficients a2,m.o(i) and C2,m.fc(i) associated with 
m ~ —2 rapidly become the dominant products. In the 
case of the orbital inclination i = 85°, they are approxi- 
mately two orders of magnitude larger in absolute value 
than the products associated with the azimuthal number 
m = 0. The products of the coefficients associated with 
m — 2 are generally two or more orders of magnitude 
smaller than the products associated with the azimuthal 
numbers m = — 2 and m = 0. 

8. Tidally induced radial-velocity variations in a 
5 Mq zero-age main sequence star 

We have applied the expressions derived in the previous 
sections to a binary consisting of a 5 Mq zero-age main 
sequence star and a 1.4 companion which we approxi- 
mate by a point mass. We considered orbital periods rang- 
ing from 2 to 6 days and the orbital eccentricities e = 0.3 
and e = 0.5. The rotation of the star is assumed to be 
synchronised with the orbital motion of the companion in 
the periastron of its relative orbit. In addition, we have 
set both the longitude of the periastron to and the time of 
periastron passage r equal to zero. For the limb-darkening 
coefiicient we adopted the value u = 0.36, which is appro- 
priate for main-sequence B-type stars. 

The observed amplitude of the tidally induced radial- 
velocity variations is displayed in Figs. |^ and 4 as a func- 
tion of the orbital period. The orbital inclination varies 
from i = 0° in the upper panels of the figures to i = 90° 
in the lower panels. The numerous peaks appearing at 
shorter orbital periods correspond to resonances of dy- 
namic tides with free oscillation modes. Bearing in mind 
that, in the perturbation theory used in Paper I, the rel- 
ative frequency difference s is assumed to be of the order 



B. Willems and C. Aerts; Tidally induced radial- velocity variations in close binaries 



7 



2 p 

:. 

-2 - 

1 

-4 - 
-6 



i = 5°. e = 0.3 



2 - 

- 

-2 

-4 - 

-6 - 



1 = 85°, e = 0.3 
* ■ . • . 

• • . 



M 
O 



-10 - 
-12 - 
-14 2^ 




m = 
m = 2 



10 
k 



15 



20 



M 
O 



-10 - 
-12 - 
-14 2^ 




m = -2 
m = 
m = 2 



10 
k 



15 



20 





2 


1 1 1 1 1 


1 1 1 1 1 1 1 1 1 1 
1 = 5°, e = 0.5 




2 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 
i = 85", e = 0.5 


CD 
M 

e 

w 

a 



-2 
-4 


«... 

■ 


■ ^ 

. . ; : 


E 

o 




-2 
-4 


■ . : 

• ■ . . • • • . . 

■ . • • . 


o 

e 
a 


-6 


* * * ^ 




o 
E 

ca 


-6 




o 








c 






log. 


-8 


- m = -2 : • 
m = : ■ 




o 


-8 


- m = -2 : • 
m = : ■ 




-10 


m = 2 : • 

III 


1 1 1 1 1 1 1 1 1 1 1 




-10 


_ m = 2 : ■> 




5 


10 15 20 




5 10 15 20 



Fig. 2. Logarithmic representation of the absolute value of the products a2^m,oii) C2,m,fc(e) for the orbital inclinations 
1 = 5° (left-hand panels) and i = 85° (right-hand panels), and the orbital eccentricities e = 0.3 (top panels) and 
e = 0.5 (bottom panels). 



of £t, the calculations near the resonances are restricted 
to values of e larger than or equal to 0.1 in absolute 
value. 

When the orbital plane is observed face-on, i.e. when 
1 = 0°, the tidally induced radial- velocity variations are 
determined solely by the dynamic tides associated with 
the azimuthal number m = 0. For higher values of the 
orbital inclination, additional contributions arise from the 
dynamic tides associated with the azimuthal numbers m = 
±2, although the contributions associated with m = 2 are 
usually negligible. The amplitude of the tidally induced 
radial- velocity variations is largest when the orbital plane 
is seen edge-on, i.e. when i = 90°. The resonances are then 
almost exclusively due to the dynamic tides associated 
with the azimuthal number m = —2. 

For close resonances, the amplitudes of the tidally in- 
duced radial- velocity variations can be quite large so that 
they are certainly detectable in observations. They in- 
crease with increasing values of the orbital eccentricity 
due to the smaller periastron distances associated with 
higher orbital eccentricities. The number of resonances is 



also larger for the orbital eccentricity e = 0.5 than for the 
orbital eccentricity e = 0.3 due to the larger number of 
terms that has to be taken into account in Expansion ^ 
of the tide-generating potential. For some cases our results 
point towards very high amplitudes which are above the 
typical sound speed in the atmosphere of the star. This 
would imply shock waves to occur, but our linear theory 
breaks down at these high velocities. The study of the non- 
linear effects is currently beyond the scope of our investi- 
gation and the applications are limited to radial-velocity 
variations that can be well described by the linear theory. 

A detailed representation of the effects of resonances 
in the range of orbital periods from 3.3 to 3.7 days is given 
in Fig. H for the orbital eccentricity e = 0.5 and the orbital 
inclinations i = 20° and i = 70°. The strongest resonances 
observed give rise to radial- velocity variations with ampli- 
tudes well above 5kms~^. Several orbital periods lead to 
resonances with two oscillation modes simultaneously. The 
peak in the amplitude near Torb = 3.37 days, for instance, 
is caused by simultaneous resonances with the oscillation 
modes ^-^d g^^. 
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In Fig. ^, the tidally induced radial- velocity variations 
corresponding to the points labeled (a) - (e) are shown as 
a function of the orbital phase. The central panel shows 
that for an exact resonance - point (c), the radial- velocity 
variations are independent of the orbital phase. The rea- 
son is that the tidally induced radial-velocity variations 
given by Expansion (^) are dominated by the term cor- 
responding to the resonant dynamic tide. Outside the res- 
onance - cases (a) and (e), the radial- velocity variations 
are markedly less regular functions of the orbital phase 
since they are now the result of a superposition of purely 
non-resonant dynamic tides. 

9. Comparison with observations: HD 177863 

The bright slowly pulsating B star (SPB) HD 177863 



was d iscovere d to be a radial-velocity varia ble by Eggen 
(1977|) , while [Waelkens fc Rufcner (1985|) detected in- 



The observed amplitude of the tidally induced 
ilocity variations as a function of the orbital pe- 
the 5Mq ZAMS stellar model and the orbital 
uty e = 0.5. 



trinsic photometric variability in this object. Later on 
Waelkens (199l| ) classified the star as an SPB. De C at 



et al. (2000) showed this star to be a single-lined spec- 
troscopic binary and determined the orbital parameters 
of the system. They found an orbital period Tdb = 
11.9154 ± 0.0009 days, an eccentricity e = 0.60, a longi- 
tude of the periastron zu = 182°, and a time of periastron 
passage t = 2450155.78 JD. The large eccentricity implies 
that the system is either still very young or that it has un- 
dergone some episodes of mass transfer. We find no indi- 
cation of the occurrence of such episodes in the literature, 
while the relatively young age is confirmed by the study 
of |Westin (1985| ), who finds an age of 64 million years. In 
what follows, we therefore assume that the system consists 
of two main-sequence stars. From the spectroscopic anal- 
ysis and the mass function, one can furthermore derive 
that the orbital inclination is between 35 and 90 degrees 
and that the companion has a mass smaller than 2 Mq. We 
also note that De Cat & Aerts (in preparation) have found 
the primary to rotate supersynchronously with a projected 
equatorial rotation velocity between 45 and 60kms~"'^. 
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0.84059c.d~i and O.lOlOSc.d-^ The first one of these dif- 
fers less than 0.001 c.d~^ from 10 times the orbital fre- 
quency. This observation, together with the binary con- 
figuration, led De Cat (2001) to conclude that one may be 
dealing with a resonantly excited mode. 

In order to investigate if the suggestion by De Cat 



Fig. 5. Detailed representation of the observed ampli- 
tude of the tidally induced radial-velocity variations of 
the 5 Mq ZAMS stellar model for the orbital eccentricity 
e — 0.5, and the orbital inclinations i ~ 20° (dotted line) 
and i ~ 70° (solid line). 



_ 2 

> 1 

" 

> 
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(200l| ) is supported by the theory developed in our paper, 
we determined the tidal response of a 3.5 Mq stellar model 
for a fixed orbital period of 11.9154 days as a function of 
the rotational angular velocity O. The model has an age 
of 64 million years, corresponding to a central hydrogen 
abundance Xc = 0.6 and a radius Ri — 2.48 i?©. The 
limb-darkening coefficient is assumed to take the value 
u = 0.36. 

The resulting amplitudes of the tidally induced radial- 
velocity variations are shown in Fig. |^ for the orbital in- 
clinations i — 90°, i = 55°, and i = 35°. Since the expres- 
sions established in our theory are derived in the linear 
approximation, we restrict ourselves to the presentation 
of radial-velocity variations with amplitudes smaller than 
the sound speed in the atmosphere of the star. The range 
of rotational angular velocities considered for each inclina- 
tion is determined from the range of projected equatorial 
rotation velocities. The companion mass for each inclina- 
tion results from the mass function of the binary. 

A large number of resonances is found for all three 
orbital inclinations, among which a whole range are caused 
by partial dynamic tides associated with the Fourier index 
= 10 in Expansion (^ of the tide-generating potential. 
The azimuthal number of the resonant dynamic tide takes 
the value m = —2 in all these cases, which is compatible 
with a preliminary mode identification performed by De 



Cat (2001). While the photometric data did not allow him 



0.2 0.4 0.6 0.8 1 
phase 



Fig. 6. Observed tidally induced radial-velocity variations 
of the 5 Mq ZAMS stellar model for orbital periods near 
3.43 days, the orbital eccentricity e = 0.5, and the orbital 
inclinations i = 20° (dashed line) and i = 70° (solid fine). 



De Cat (2001) also performed a detailed frequency 
study on both extensive spectroscopic and photometric 
time series and found two intrinsic frequencies for the star: 



to identify the mode, he found convincing evidence from 
the variations of th e moments of the line profiles (for a 
definition see, e.g., Aerts ct al. 1992 ) to conclude that 
the observed mode is most likely a sectoral £ — 2 mode. 
We note that any resonances caused by partial dynamic 
tides associated with the azimuthal number m = would 
have led to large amplitudes of the tidally induced radial 
velocity variations for all values of the rotational angular 
velocity. 

In the case of the orbital inclination i — 35°, we in- 
dicated the positions of the resonances associated with 
m = — 2 and fc = 10 by dotted vertical lines in Fig. ^ 
We find a possible resonant excitation of sectoral £ = 2 
modes ranging from the mode (/^ up to the mode g^. For 
the orbital inclinations i — 55° and i = 90°, the radial 
orders of the modes excited by the partial dynamic tide 
associated with m — —2 and fc = 10 are even higher. 

Upon close inspection, the pattern of resonances in 
Fig. 1^ can be seen to repeat itself at equidistantly spaced 
values of the rotational angular velocity fi. The reason is 
easily identified when the forcing angular frequency ctt is 
rewritten as 



kn + mil ~ (k + l)n + m (fl ) . 

V m/ 



(37) 
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Fig. 7. Observed amplitude of the tidally induced radial-velocity variations in HD 177863, for the orbital inclinations 
i — 90° (left), i = 55° (middle), and i = 35° (right). The dashed lines in the left- and right-hand panels indicate the 
rotational angular velocities corresponding to 10% and 20% of the critical break-up angular velocity ilc- The dotted 
lines in the right-hand panel indicate resonances caused by the partial dynamic tide associated with m = — 2 and 
k = 10. 



It follows that if the partial dynamic tide associated with 
the azimuthal number m and the Fourier index k is reso- 
nant at the rotational angular velocity il, the partial dy- 
namic tide associated with the azimuthal number m and 
the Fourier index k + I will be resonant at the rotational 
angular velocity il — n/m. 

We conclude that due to the large uncertainty in the 
rotational angular velocity and the orbital inclination of 
HD 177863, a definitive mode identification is not yet pos- 
sible at this time. In addition, the amplitude of the tidally 
induced radial- velocity variations can not be determined 
in the isentropic approximation so that the inclusion of 
nonadiabatic effects in our treatment becomes desirable. 
For higher rotation rates, more resonances can also be ex- 
pected due to the rotational splitting of the stellar eigen- 
frequencies by the Coriolis force. We plan to explore the in- 
fluence of the nonadiabatic effects and the Coriolis force in 
more detail in subsequent investigations. In particular, we 
will investigate the behaviour of the predicted observable 
quantities (photometric variations and line-profile varia- 
tions) resulting from the candidate resonant modes that 
we have found in this work with the goal to confront them 
with the data. Such an iterative procedure can perhaps 
allow us to come to a definite mode identification and de- 
tailed modelling with the goal to constrain the internal 
structure parameters of the primary. 

10. Concluding remarks 

In this paper, we derived semi-analytical expressions for 
the tidally induced radial-velocity variations in a uni- 
formly rotating component of a close binary. We neglect 
the effects of the Coriolis force and the centrifugal force 
and treat both the free and the forced oscillations of the 
component as linear, isentropic perturbations of a spheri- 
cally symmetric star. We take into account the possibility 
of resonances between dynamic tides and free oscillation 



modes by means of the perturbation theory developed by 
^meyers et al. 



The amplitudes of the tidally induced radial-velocity 
variations depend on the orbital eccentricity and on the 
inclination of the orbital plane with respect to the plane 
perpendicular to the line of sight. They increase with in- 
creasing values of the orbital eccentricity and are largest 
when the orbital plane is seen edge-on. The amplitude 
is most sensitive to the values of the orbital inclination 
when 20° i ^ 70°. A similar conclusion was reached by 



Kumar et al. (1995) in an investigation on the stellar lumi- 



nosity variations associated with tidally excited oscillation 
modes. 

From the application to a 5 Mq zero-age main sequence 
star, it follows that the amplitude of the tidally induced 
radial-velocity variations seen by an observer is small, ex- 
cept when resonances occur between dynamic tides and 
free oscillations modes. The resonances enhance the tidal 
motions of the mass elements and lead to radial-velocity 
variations with amplitudes that are up to an order of 
magnitude larger than those observed outside resonances. 
Some orbital periods are seen to give rise to simultaneous 
resonances with two oscillation modes. 

The shape of the tidally induced radial- velocity curves 
changes markedly with the proximity of a dynamic tide 
to a resonance with a free oscillation mode. For close 
resonances, the tidally induced radial-velocity variations 
are almost exclusively caused by the tide involved in the 
resonance and therefore exhibit a very sinusoidal-like be- 
haviour. Outside resonance, the radial- velocity variations 
are less regular functions of the orbital phase. 

For conclusion, we applied our results to the slowly pul- 
sating B star HD 177863 and showed the possible resonant 
excitation of a high-order second-degree sectoral g^-mode 
in this star. 
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The functions f^.e{u) and fri.iiu) defined by Eqs. ( ^31 ) and 
( p^ render the influence of the adopted limb-darkening 
law on the contributions of the radial and the transverse 
component of the tidal displacement to the tidally in- 
duced radial-velocity variations. In the case of the limb- 
darkening law given by Eq. (|2l|), the functions f^,t{u) and 
frii{u) take the form 
6 

fijXu) = - — - [(1 - M)7i',2 + U7f,3] , (A.l) 



3 — M 
6 



£ [(1 - u) (7£,2 - Je-i,i) 



3 — u 

where we have introduced the abbreviation 
7i,n = / Pd^) x"- dx. 



(A.2) 



(A.3) 



For non-negativ e val ues of n, the integral in the right-hand 
member of Eq. ( A.3 ) yields 

7ri/2r(7H- 1) 



(A.4) 



2«+ir[(n-£ + 2)/2]r [(n + £ 3)/2] ' 
[lAbramowitz & Stegun|[l965i Eq. (8.14.15)]. 

In the particular case where ^ = 2, the functions f^,2{u) 
and frf.2{u) take the form 
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(A.5) 
(A.6) 



